Generating correlated (2+l)-photon in an active Raman gain medium 
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A scheme of generating controllable (2+1) photons in a double-A atomic system based on active- 
Raman-gain is presented in this paper. Such (2+1) photons can be a potential candidate to generate 
a correlated photon pair as one photon of 2 photons acts as a trigger. Our proposal is an alter- 
native approach to generate photon pairs where the frequencies are tunable for subsequent atomic 
system experiments compared to spontaneous parametric down-conversion case. Compared to other 
schemes of generating correlated photon pairs, our scheme exhibits several features due to the exploit 
of the stimulated Raman process and injection-seeding mechanism. 
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I. INTRODUCTION 

Spontaneous parametric down-conversion (SPDC) is a 
widely used method for producing correlated and en- 
tangled photon pairs. Due to the spontaneous emis- 
sion nature of the process, SPDC-based photon sources 
usually have limited applications because of broadband 
spectrum, low efficiency, short coherence time and co- 
herence length. Different approaches to generation of 
photon pairs have been demonstrated experimentally by 
a few groups 0, i, S i, 1, 1, 0, 0. A few years 
ago, Kimble and co-workers working with a magneto- 
optical trap fl[ and Lukin and co-workers working with 
hot atoms |2j have shown the generation of nonclassi- 
cal photon pairs. Matsukevich and Kuzmich 0] real- 
ized the nonclassical photon pairs using the two dis- 
tinct pencil-shaped components of an atomic ensem- 
ble. Harris and co-workers [J demonstrated the gener- 
ation of counter-propagating paired photons in a double- 
A four wave mixing (FWM) scheme where the anti- 
Stokes/FWM field is generated under the electromag- 
netically induced transparency (EIT) condition. With a 
modified scheme Kolchin et al. @ showed the generation 
of narrowband photon pairs using a standing wave pro- 
duced by a single intense driving laser. Thompson et al. 
Q reported the generation of narrowband pairs of pho- 
tons from a laser-cooled atomic ensemble inside an opti- 
cal cavity. Du et al. 0| reported the production of bipho- 
tons in a two-level atomic system. Pan and co-workers [§J 
demonstrated polarization-entangled photon pairs using 
only simple linear optical elements and single photons. 
A key element of these experiments H, H, EL IE IE 0] is 
that one photon is generated by the spontaneous Raman 
emission process whereas the second photon is generated 
via an ElT-assisted two-wave mixing process. Alterna- 
tively, the process can be viewed as a FWM generation 
with one spontaneous emission step in the wave mixing 
loop. 

In this paper, we present a far-detuned, active Raman 
gain (ARG) scheme 0, HE HH, E2] for the generation 
of a group of correlated (2+1) photons. Our goal is to 



present a different and alternative approach to generate 
photon pairs where the frequencies are tunable for sub- 
sequent atomic system experiments compared to SPDC 
case. Especially, when a perfect single source is on de- 
mand our scheme will have potential applications. At 
the first glance the scheme reported here is very similar 
to those studied in Refs. 0, 0, II 0, 1, 1, . However, it is 
operated under a very different principle. Instead of rely- 
ing on a spontaneous Raman process to generate the first 
photon, we inject a probe photon which leads to stimu- 
lated Raman generation. Two features are evident in this 
scheme: (l)injection seeding of the probe photon leads to 
highly directional generation of probe photons and FWM 
photons, (2) stimulated Raman process ensures that the 
more photons appear in the same frequency mode and 
the gain can be easily controllable. The first feature ef- 
fectively increases detection efficiency in comparison with 
a process that relies on spontaneous emission of the first 
photon. Due to stimulated emission the generated probe 
field has the band width and frequency characteristics 
exactly as the inject-seeding field. The second feature 
can lead to the generation of photon-number Fock state, 
which will be explained in Sec. [V] The precisely con- 
trollable directions of photons due to injection-seeding 
mechanism in our scheme offer a larger flexibility for dif- 
ferent applications than the undetermined directions of 
photons in SPDC case. 

Carefully choosing the detuning and intensity of the 
second pump field, our scheme generates a group of (2+1) 
photons where two photons are in the same probe fre- 
quency mode whereas one photon is in the FWM fre- 
quency mode. The schematic of the process is shown in 
Fig. [1] We consider an ensemble of identical lifetime 
broadened four-state atoms initially prepared in their 
ground states |1) by suitable optical pumping method. 
A cw pump field (-Epi) couples the ground state |1) to 
an excited state |2) with a large one-photon detuning 5\. 
A second strong cw laser field {Ep2) couples the state |3) 
to an excited state |4) with the detuning A + S4. A single- 
photon probe field E\ (with central frequency W3), which 
couples the state |2) to a lower excited state |3), is then 
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FIG. 1: (Color online) (a) Schematic atomic level diagram. 
.Epi(wi) and Ep2(u>2) are the pump fields, and Ei{ujz) is the 
probe field. E2(u>4) describes the process of another photon 
(0^4) generated by FWM. (b) Schematic of the experiment for 
2+1 photons generation. 



introduced to the medium. Consequently, a two-photon 
Raman transition can occur, causing an atom to absorb 
a pump photon and emit a photon into the probe field 
frequency mode. The presence of the second pump field 
subsequently pumps the atom to state |4) and a FWM 
process occurs, leading to the generation of a new photon 
at the frequency L04 = lo\ — 103 + L02- In such a process two 
photons at 0J3 as well as one photon at L04 are simulta- 
neously gen erated, yielding a group of (2+1) correlated 
photons [13]. 

The paper is organized as follows. In Sec. [IT] we first 
derive the equations of motion for atomic dynamics and 
field propagation equations. With analytical solutions we 
first examine a limiting case where our scheme reduces 
to a well-known single A gain scheme. We then investi- 
gate group velocities of various propagation modes and 
corresponding propagation parameters for the injection- 
seeded double- A scheme. In Sec. lIIIl we consider the case 
where a single-photon quantum probe field is injected 
into the medium and we investigate the generation of 
the two probe photons and one FWM photon under the 
condition of weak gain. We also analyze the two-photon 
intensity correlation function and coincidence count rate. 
We also describe the case where the input probe field is a 
coherent state and obtain the single-photon-added coher- 
ent state. In Sec. IIV1 we show the conversion efficiency 
of our scheme and compare the efficiencies between our 
scheme and SPDC case. In Sec. [V] we describe how to 
make our scheme be an alternative approach to generate 
the photon pair and the possible applications. In Sec. lVIl 
we discuss possible complications due to various ac Stark 
effects. We also discuss the bandwidth of the Raman 
gain and the efficiency and bandwidth of single photon 
source. A summary is given in Sec. I VIII 



II. MODEL AND ANALYSIS 
A. Theoretical model 

In the electric-dipole and rotating-wave approxima- 
tions, the interaction Hamiltonian of N identical four- 
state atoms interacting with laser fields (see Fig. 1) is 
given as 

(1) 

where Eq + ^ is the slowly varying quantum field operator 
(q = 1,2 are for probe and FWM field, respectively), 
E { - ] = {E { + ] )\ and & 3 kl = \k) ]3 {l\ (k,l = 1,2,3,4). 
In addition, 2Qi = ^2iEpi/h and 2Q2 = A*43-Ep2/^ 
are Rabi frequencies of pump fields with \i nn i being the 
dipole matrix element for the \n) — \n ) transition and 
Ep\ t p2 being the amplitudes of the two classical pump 
fields, gi and g 2 are the atom-photon coupling constants, 
and u> m and k m (m = 1, 2, 3, 4) are the carrier frequency 
and wave vector of the m — th optical field, respectively. 

We define continuum atomic operators 0^ by sum- 
ming over the individual atoms in a small volume V, and 
introduce slowly varying atomic operators 5^: G\2 = 
e- luJlt a 12 ,a 13 = e^-^aia = e^-^a^, a u = 

e-^CTu, CT 23 = e^ 3t a 23 ,cr 2 4 = e-^-^*CT 2 4 = 
e - l (^ 2 -^ 3 )t 5 - 24jCr34 = e -^ 2 t~ 34j the equations f or tnc 

atomic operators in the Heisenberg picture are 

0-21 = -(721 - ^1)0-21 + iflle~ lkl ' r (a 2 2 - oil) 

+ ig 2 Et ) e- i ^ f a 24 -ig 1 E[- ) e- i ^ f! a 31 , (2a) 

034 = -(734 + ^2)0-34 - iQ. 2 e lk% ' r {au - cr 33 ) 

- i gi E{ +) e^- p a24+ig2Ei +) e^- p a 31 , (2b) 

031 = -(731 - *A)a 3i + itlle^-%2 - in 2 e l ^- p a 41 

- ig 1 E[ +) a 21 e i ^+ig 2 E 2 - ) e- i ^a 34 , (2c) 
ct 41 = -(7 41 - 164)0-41 + i0^e~ ifcl ' F o-42 - i^l*2 e ~ %k ^ ? 

- ig2E ( -\a ll -cj AA )e-^ 1 (2d) 

032 = -(732 + ^3)032 + ifiie ifcl "*CT-3i - ^2e 4fc2,r "o 4 2 

- z3i^i +) (o22-o 33 )e^- F , (2e) 
042 = -[742 - «(^i - <5 4 )]ct42 + i^io- 4 ie 4fcl ' r - 1^2032 

x e-^ + igtE^a^e^ 

- ig2E { 2- ) ~a 1 2e-^\ (2f) 

where 5i = uj 2 i - 0J\, S 2 = w 43 - u 2 , S3 = cj 2 3 - ^3> $4, = 
LU41 — to 4 , A = S3 — 81 = 62 — 84, and u>ki — tuk — is 
the frequency of the \k) <-> \l) transition. 7^ (7^ = 7;^.) 
is the dephasing rate between state \k) and state 

We consider a pencil-shaped atomic ensemble and as- 
sume the atomic sample is optically thin in the trans- 
verse direction. Using the slowly varying envelope and 
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unfocused plane-wave approximations, we obtain the fol- 
lowing propagation equations for the quantum field op- 
erators: 



d d 



d d 



(+) 



(z,t) 



ig 1 Na, 



32 • 



(3) 



-—ik,A 



= -ig 2 Na Al , (4) 



where N is the total number of atoms in the atomic en- 
semble. Under the condition that the probe field is much 
weaker than the pump fields, and using the assumption 
that pump fields propagate without depletion, Eqs. (j2a|) 
and J2fj can be evaluated adiabatically, 



ik\ -z 



«742 



Considering = 



'2i ~ Si + i-ni' 42 ~ s 1 

(5) 

a 44 = ^34 ~ 0, an d substi- 
tuting these values of the density matrix elements into 
Eqs. (Pc)) - Pe|) . we obtain the following set of coupled 
equations: 



5-31 = -(731 - «A)ct 3 i + inie~ ikl - z d Z2 
(741 = —(741 - ^4)5-41 + ift*e~ 



CT 42 



zf2 2 e 



&31 - 252^2 [z,t)e , 



, CT32 = -(732 + ££3)032 + iSli<7 3 ie lfel ' 



(6) 

Substituting Eq. (J5J) into Eq. ©, we make the Fourier 
transform and obtain the solutions for cr 32 and 041 in the 
frequency domain 



^32 



E41 — 



diD{uj 

ff2^1^2 

g 2 [(w + tfc)(w 



ik% -z 



(u> + d 5 ) 4(z, _ w ) e *&+£»-£0-* (7) 



2 {tu-d 3 ) ei(z,w)e~ l(fcl+ ' £2 ~' £3) - ? 



D(u) 



4(z,-u>)e- ik ^, 
(8) 



where di = Si + 1721, d 2 = A + 1734, d 3 = 5 3 - Z7 32 , 
c?4 = £4 + 1741, d 5 = Si + 4742, and D(ui) — (u + 
d 2 ){d 3 - uj)(lu + di) + |Sli| 2 (a; + d A ) + \Q 2 \ 2 (lu - d 3 ) + 
\riiil 2 \ 2 (uj + 2dg)/<i|. Here we ignore some small terms 
under the condition <5i(<5 3 ) 3> S4, A,£li,£l 2 . In addi- 
tion, e 3 {z,uj) = f^E^iz^ty^dt, and E 32 (z,u) and 
£44(2,0;) are Fourier transforms of a 32 (z,t) and 044(2, f), 
with to being the transform variable, respectively. Mak- 
ing the Fourier transform of Eqs. ^ and Q and using 



Eqs. and © we obtain 

d iuj . . A . 

— ei(z,o;) = — £4(2, u>) + zDi(a;)ei(z, w) 
az c 

+iZ^)402,-u;)e iA ^, 
^-4(*»-w) = — 4(«>~w) +i£' 3 (w)ei(z,c i ;)e- 4Afc 



+iD 4 (w)4( 2; I -^) 1 



(9) 

•Z 

(10) 



where 



didsUyuJ) 

/f 12 nin 2 (w + «z 5 ) n , x Jfi2i2!n5(w-d 3 ) 

^M^J — „/ \ 1 5-^3^] = - 



D{uj)d 5 ' J ' 7 D(cj)di 
K 2 [(u> + d 2 )(u - d 3 ) - l^l 2 } 
D{u) 



(11) 



in which K x = N\ 9l \ 2 /c, K 2 = N\g 2 \ 2 /c,Ki 2 = Ngig 2 /c, 
and phase mismatch Afc = fc4 — fc 3 +fc 2 — fci. Equations ([9]) 
and (JTU]) can be solved analytically for arbitrary initial 
conditions. For simplicity, we let the phase mismatch 
Afc = 0. The solutions of Eqs. © and (JTUJ) are as follows 



ei(z,u>) 



el(z,-uj) 



U+ - U. 



■[{U+e lX+z - C/_e iA - z )e!(0,w) 



+ [/ + t/_(e 4A - z -e 4A+2 )4(0,-w)], (12) 
1 [(e lA + z -e jA - z )ei(0,w) 



U+ - U- 



(U+e 



iX- z 



[/_e iA + z )4(0,-c)], (13) 



where 



A± = ~ + hD 1 (u) + Di(w)TD s (w) 
c I 



2D 2 (uj) 



Di(u) - Di(cj) T D 5 (uj) 



D 5 (lj) = >/(Di(w) - Di(u)) 2 +4D 2 (w) J D 3 (w).(14) 

We focus our attention on the adiabatic regime [3], 
where A± and U± can be expanded into a rapidly con- 
verging power series of dimensionless transform variable 
£ = LUTp where r p is the pulse duration of the weak 
probe field. In this regime [HI, U± = W± + 0(C) and 
A± = (A±)^=o + Q/T p V g ± +C(C 2 ) can accurately describe 
the FWM generation and propagation process. The ex- 
pansions can be well valid close to the central frequency 
component of the field. The inverse Fourier transform of 
Eqs. dHJ) and (13]) is given by 



E[ +) (z,t) = [A 1 E^(n+)-A 2 S^'( n +)] 

+ [A 2 E { -\^)-A 3 E { +\^)y- z llb) 
E { -\z,t) = [AE[ + \ V+ )-A 3 E ( -\ V+ )}e^ 

+ [AiE ( 2 -\r 1 ^-AE{ +) (r 1 ^- z , (16) 
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FIG. 2: (Color online) The schematic diagram of fields E\ and 
E2 propagation. Two propagation modes are characterized by 
the different group velocities V g± where one mode gains and 
the other mode attenuates. 



where A = 1/(W+ - W-), A x = W+A, A 2 = W+W-A, 
A 3 = W~A 7 and r)±=t- z/V g ±, f3± = i(A±)w=o and 



W± 

At 
1 

vT 



2D 2 (0)/(D 4 (0)-£>i(0)t£> 5 (0)), 
i[£>i(0)+X> 4 (0)T ^5(0)1/2, 



9± 



D lp (0) + D 4p (0) T D 5p (0) 



(17) 
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D p (0) 
D lp (0) 



I no 



KM 2 



d 3 d 4 - 
{d 4 d 5 



d?d 4 



ir> 2 



d 2 d 3 , 

-Nil 2 ; 



D(0) 1 



did 5 D(Q) 

D 2p (0) = -ifi 2 fiin 2 (£)(0) + d 5 ^ P (0))/d 5 ^ 2 (0), 

^ 12 ^f>*£ p (0)OD(0) + d 3 ^ P (0))/di^ 2 (0) 

#2 



3p(0) 
#4p(0) 



^5 P (0) 



/}J(0) [(da - rf 3 )£»(o) + ^ P (0)(d 2 d 3 + l^l 2 )], 
[(£>i(0) - D 4 (0))(Ui p (0) - D 4p (Q)) 



-2D 2p (0)D 3 (0) + 2D 2 (0)£> 3p (0)]/ A /£> 5 (0). 

(18) 



Here ^(0) = Dy^)]^, (j 



1. 



,5) and D(0) 



-D(w)| w =o- Equations (fT5|) and (|16| indicate that in 
general each frequency component of the probe and the 
generated FWM fields contains two propagation modes 
(wave packets) that travel with different yet individually 
matched group velocities. In addition, both propagation 
modes (wave packets) retain a pulse shape identical to 
that of the input probe field in the adiabatic regime [lj| • 



B. Analysis of group velocities and propagation 
parameters 

Before a detailed analysis we first consider a limiting 
case where J7 2 = and g 2 — 0. In this limit, our scheme 
reduces to that of a single A scheme. Using Eq. ^ we 
immediately obtain [note that when g 2 = 0, D 2 (lo) = 0], 
thus ii(z,u>) = ei(0,uj)e~ lX ° z where Ao = ui/c + -Di(w). 
This leads to a single group velocity 
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FIG. 3: (Color online) Group velocities Re[V s+ /c] (solid line) 
and Re[V 9 _/c] (dashed line) as function of dimensionless two- 
photon detuning Ar p for K\ — K2 = K12 = 1 x 10 9 /(m-s), 

[Oil = 7, |n 2 j = 57, t p = 10 /is, 61 = 1007, 64 = 0.17, 

731 = 3 x 10~ 5 7, 721 = 723 = 741 = 743 = 7 = 10 MHz, and 
742 = 27. The inset shows a magnification of the region for 
small V g+ and V g _ . 



Alternatively, one can also obtain this result using Eq. 
(12). Note that with Q 2 = and g 2 = 0, we have 
D 2 {uo) = D 3 (w) = D 4 {uj) = 0, D 5 {uj) = 
U+ = U- = 0, and A_ = A = uj/c + D\{uj) and 
e\{z,w) = €i(0, w)e~ lX ° z . Equation (fT^|) is the exact 
same superluminal group velocity obtained by Payne and 
Deng [9( in an atomic amplitude treatment of a single A 
active- Raman-gain scheme. 

In general, for an injection-seeded double-A active gain 
system presented here, each field consists, as indicated 
in Eqs. (fl"2"|) and (fill)) , of two propagation modes char- 
acterized by the eigenvalues A±, and therefore different 
group velocities V g± and decay rates. In the adiabatic 
regime, the propagation of the fields is simply governed 
by Eqs. fT5]) and (JTHJ). The schematic diagram of fields 
Ei and E 2 propagation is showed in Fig. [21 For a short 
propagation distance where the decay of individual mode 
is not significant, these modes will overlap each other. 
Theoretically, for a sufficiently long propagation distance, 
however, modes with different propagation velocities will 
separate and the fast decay components decay out com- 
pletely. Consequently, one obtains a pair of well- matched 
waves traveling with the identical group velocity and have 
the identical decay behavior. 

In Fig. [3] we plot group velocities of each propagation 
mode (wave packet) as a function of dimensionless two- 
photon detuning At p (i.e., the ratio of two-photon detun- 
ing to the bandwidth of the probe field) for a typical cold 
alkali vapor where states |2) and |4) belong to the same 
hyperfine manifold. When the two-photon detuning is 
large, i.e., Ar p > 22.4 or Ar p < —24.3 for this specific 
example, one wave-packet mode travels with a negative 
(superluminal) group velocity and the other wave-packet 
mode travels with a positive (subluminal) group veloc- 
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FIG. 4: (Color online) The coefficients Re[/3±] versus the di- 
mensionless two-photon detuning At p . The parameters are 
chosen as those using in Fig. [3] 



ity that can be substantially smaller than the speed of 
light in vacuum. In this case, the fast and slow modes 
separate quickly even before the differential decays be- 
come significant. Consequently, one obtains two group 
velocity-matched probe-FWM field pairs arriving at the 
detector at a delayed time [15]. When the two-photon 
detuning is small both propagation modes travel with 
positive group velocities. Under the condition where the 
two group velocities are equal, one may obtain only one 
pair of group- velocity matched fields. 

To further understand the propagation behavior of 
probe and generated fields, we examine the propaga- 
tion parameters (3± in Eqs. (fT5|) and (flU)) . The case of 
Re[/3±] > corresponds to gain whereas Re[/3±] < in- 
dicates field attenuation through propagation. In Fig. [4] 
we show these propagation parameters as functions of 
the dimcnsionless two-photon detuning Ar p . In the re- 
gion (—10.8 < At p < 8.8) where both propagation modes 
(wave packets) travel with positive group velocities one 
propagation mode experiences amplification whereas the 
other propagation mode is attenuated. Note also that 
near Ar p w 0.305 the signs of propagation constant 
change and the amplified (attenuated) mode becomes at- 
tenuated (amplified). These features are very different 
from the conventional EIT based FWM schemes where 
probe attenuation and FWM gain always occur. This is 
precisely due to the fact that EIT process is based on the 
weak absorption of the probe field and stimulated genera- 
tion of the FWM field. In the case of active-Raman-gain 
medium, the probe field serves as an injection seeding 
source and it works in a stimulated emission mode. Con- 
sequently, simultaneous gain to both probe and FWM 
fields is possible. Indeed, if parameters are chosen such as 
Re[/?+] > and Re[/3_] < 0, then after a sufficient prop- 
agation distance, both the probe and FWM fields have 
the same gain feature characterized by Re[/3+] > 0. The 
energy that supports this bi-field increase comes from the 



two CW classical fields Epi and Ep2- 

III. GENERATION OF GROUPED AND 
PAIRED PHOTONS 

In this section, we consider the case where the quantum 
state of the injected probe field corresponds to a single- 
photon wave-packet state [H, U3] : 

/oo 
dwP p i(n7 + w')fit(<i/)|0}, (20) 
-oo 

where the amplitudes P p \ (w + to ) are normalized such 
that du \P p i(zu + uj )| 2 = 1 and the w is central 
frequency of wave packet. Hence the initial state for the 
system is 

IVO = |1> 3 10> 4 , (21) 

where subscripts 3 and 4 denote single photon wave 
packet of central frequency uj 3 and u>4, respectively. In 
general, the generated quantum state of the system at 
time t can be expanded in terms of boson Fock space as 

IVw) = ^a nm (t)\n} 3 \m) 4 , (22) 
n m 

where n and m denote the photon numbers of the fields 
with the central frequency tt>3 and ui±, respectively and 

E„mKm(*)| 2 = L 

When the generated photon numbers {n, m} are small 

enough in the low gain case, we can work out the coeffi- 
cients {a nm } in terms of the moments of the probe and 
FWM field operators with the forms 

E^ ) (0),---)\^ out ) = ^ m \F(Et ) (L), 
Ei + \L),E^\L),E^(L),...)\i> in ), (23) 

where F(- ■ ■ ) denotes the combinations of products of 
the field operators required to calculate the coefficients 
{ctnm}- (0) an d E^^L) are the quantum field oper- 
ators at the entrance z — and at the output end z = L, 
respectively. 

A. Generation and analysis of (2+1) photon group 

To generate a (2+1) -group of photons we consider 
the weak gain limit in which the injected single-probe 
photon will generate only one photon in the probe fre- 
quency mode by stimulated Raman process, accompa- 
nied by another photon in the FWM mode. In this case, 
we only need take account of the low photon numbers: 
n, m € {0, 1, 2} in the expansion of Eq. (j2"2")l . As a result, 
the final state can be written as 

= a M |1> 3 |0) 4 + a 20 |2> 3 |0) 4 + a 21 |2) 3 |1) 4 . (24) 
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FIG. 5: (Color online) The ratio of ln|a32 |/ln|a2i | versus the 
z/z m for At p = — 1 and z m = 1 cm (solid line). The other 
parameters are chosen as those using in Fig. [3] The dashed 
line is for the case the ratio is 2 (SPDC case). 



Here, the physical meaning of each term is very clear, 
aio represents the probability amplitude of the injected 
single-probe photon not initiating a stimulated emission. 
a20 describes one photon generated through the stimu- 
lated emission in the probe mode but without photon 
generated in the FWM mode. a 2 \ is the probability am- 
plitude of one photon generated through the stimulated 
emission in the probe mode, and simultaneously one pho- 
ton generated in the FWM mode. The second term in 
Eq. (|2~4")1 only exists for a pump field fl 2 which is too weak 
to excite the FWM process. When the second pump 
field is strong enough to generate the FWM photon , 
this amplitude a 2 o tends to zero, and the second term in 
Eq. (j2"4")) vanishes. Further we assume this case to ensure 
the generation of a (2+l)-group of photons. Then the 
state vector simply reduces to the form 



m)) = |aio||l>3|0>4 + e*|aai||2> 3 |l>4. 
With the help of Eq. ([2*5]) , one can work out 



(25) 



aio 



a 2 i 



4 \A\*\eP+ L 
= 2|A|V+ 
= arg[(Ai 



.4 



0-LI2 



3& 

./3-i|2l p2 



\A 2 



pl (t-L/V g -(z-L)/c), 



eP- L \ 2 P*(t-L/V g -(z-L)/c), 



pi 

f3-L\ ( B + L 



A 3 e^)(e 



3 p - L )*A*], 



(26) 



where we consider the case with the group velocities 
V g+ = Vg- = Vg, and the pulse is assumed to main- 
tain its shape with the profile P p i(t — L/V g — (z — L)/c) 
during propagation. It is worth of pointing out that for 
a slightly higher gain, the higher order terms, such as 
a 32 1 3) 3 1 2) 4, describing the multiphoton processes will ap- 
pear in Eq. (|25|) . Similarly, multiphoton processes also 
exist in the SPDC case in the regime of high gain. For an 
ideal SPDC case, the quantum state of photons can be ex- 
panded in the form |^>) = |0) + g\U, l s ) + g 2 \2 i} 2 S ) H , 

and the probability amplitude of two-photon pairs re- 
duces in the square law for g << 1. For comparison, we 
numerically estimate the corresponding probability am- 
plitude ICK32 1 in our case. The result shown in Fig. \S\ in- 
dicates that the probability amplitude |«32| has a faster 



reduction than the square law for a length z > 2.5z m . 
This means that our scheme using stimulated Raman 
process with injection-seeding mechanism is slightly bet- 
ter than the SPDC case for compressing the multiphoton 
processes, and hence appropriate for generation of single- 
photon pair. 



B. Two-photon intensity correlation function and 
coincidence count rate 

To show the time correlation properties of the gener- 
ated photon pairs, we now work out the Glauber intensity 
correlation function between the paired photons (uj 3 , 104) 
with a time delay r^, 



= (Er'mf'it + r^'it + r^Erit)). 

(27) 

Note that the correlation function is calculated for the 
state of two photons (one in probe mode and the 
other in FWM mode), not the state of three pho- 
tons. The values ( ) is average of the initial state 
\i>in) = |1) 3 |0) 4 = |0) 4 J duP pl (-u)e\(0,-Lo)\0) 3 where 
Ppi(w) = J^dte^Ppiit) is the pulse shape 

function of the input photon. 

Using the fields operators of Eqs. (fT2|) and (fl"3]) . one 

( 2) 

can work out G ' E _ E (Td). In order to denote the result 
simply, Eqs. I|12[) and (JT3J) can be written as 

h{z,u) = iJiei (0 J o;) + 5i4(0,-w), (28) 
4(z,-w) = i? 2 e 10 (0,^) + S 2 4o(0,-^), (29) 

( 2) 

and the intensity correlation function G Ei _ E2 (jd) has 
been derived as 



•A+), 



G 



dhj\ dio 2 dio 3 dtUi e 



— iuit 



E 1 -E 2 ( Td ) ~ 
xe -» 3 (i+r,!) e -i U 3((+Ti) e -™4^t( 2i _ Wl ) 

Xel(z, -uj 2 )£ 2 (z,U3)ei(z,u>4)) 
= [J d W P p 2 1 M|i? 1 | 2 + J dLo\Si\ 2 ] J du\R 2 \ 2 

+ J duP 2 AR 2 \ 2 J du\S 1 \ 2 + \ j due WTd SlS 2 \ 2 

+ J due-^P^i^R^ x J dcue-^SiS; 

+ J due^P^MRiRZ x J dcje WTd S* 1 S 2 

« Gg(0)Gg(0) + | J dcje i ^S* 1 S 2 \ 2 

+ J duje- luJTd P 2 1 (uj)R* 1 R 2 x J duje- lUTd SiS* 2 

+ J due^P^i^RiRi x J due luJTd S* 1 S 2 . (30) 

The second-order normalized intensity correlation func- 
tion is g ( S_ E2 (r d ) = G^l E2 (T d )/G { ^(0)G^ 2 (0) where 



FIG. 7: (Color online) The schematic of generating single- 
photon added coherent state (SPACS). 
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FIG. 6: Coincidence count rate in a 1 ns bin versus the delay 
time Td for Ar p = —1 and z/ z m — 5 at (a) K — 2 x 10 8 /(m-s) 
and (b) K — 3 x 10 9 /(m-s). The other parameters are chosen 
as those using in Fig. 



consider injection of a weak coherent state |a}3 into the 
medium instead of a probe single photon 1 1) 3 . Here, we 
also consider the weak gain limit where only one pho- 
ton will be generated into the probe field by stimulated 
Raman process, accompanied by another photon in the 
FWM mode. With such an arrangement, we have the 
output state 

M =/3i|a>3|0>4 + /3 2 M> 3 |l>4. (31) 

When only a single photon in the frequency W4 is de- 
tected, the state \if>) will collapse to the SPACS \a, 1). 
The schematic is shown in Fig. [71 Again, our scheme 
using stimulated Raman process in atomic ensemble of- 
fers a direction-controllable way to generate a SPACS, 
compared to the present method demonstrated in the 
nonlinear crystal [20(. 

IV. CONVERSION EFFICIENCY 



the peak value of normalized g E E 3> 1, the antibunch- 
ing nature, which implies a finite time delay for the emis- 
sion of the second quantum field photon (W4). In other 
words, a nonclassical photon pair composed of frequen- 
cies L03 and W4 is generated in atomic vapors. 

The coincidence count rate R c (Td) is obtained 
from the intensity correlation function as R c (Td) = 
e 2 ATGg'_ & (Td) with a bin size AT = 1 ns much smaller 
than the correlation time. The factor e accounts for 
the photon counter efficiency, the filter transmission, and 
fiber coupling. In Fig. [6l we show the coincidence count 
rate R c (Td) in a 1 ns bin depending on the K. When 
K is small, which occurs, for example, at low optical 
depth, the atomic system behaves like a single atom. In 
such a regime the intensity correlation function reveals 
the damped Rabi oscillations. By increasing the opti- 
cal depth of the atomic sample, it becomes possible to 
achieve the E1-E2 correlation function with shorter time. 

C. Single-photon added coherent state 

Another possible extension of the stimulated Raman 
process with injected seed is to generate a single-photon 
added coherent state (SPACS) [19| . For this purpose, we 



In this section, we come to evaluate the output conver- 
sion efficiency of our photon pair. The intensity of the 
fields is given by [2l[ 

I = 2n(s /no) 1/2 \E\ 2 . (32) 

where e = 8.85 x 1CT 12 F/m, /x = 4tt x 1CT 7 H/m, 
and n is the refractive index, and E is measured in V/m. 
First, we calculate the intensities of pump fields Ep\ and 
Ep2- The relations of Rabi frequencies of pump fields 
with the transition matrix element are 2S1 1 = fi2iEpi/h 
and 2f^2 = [i^Ep^/h, so the intensity of Ipj in the un- 
depleted pump approximation (Ipj(z) = lpj(0)) is 

I Pj (0) = 8n^£o/Vo) 1/2 |MV^| 2 , 3 = {1, 2}, (33) 

where /ii and /i2 denote /U21 and /X43, respectively. 

Next, the intensities of generation fields E\ and E2 at 
the boundary L are 

I Ej =I 0j (E^E^), j = {1,2}, (34) 

where I 0j = ( v //io;2+i/2£o^cffCTp) 2 (j = 1,2) [13] and 
A e ff is the effective beam cross section, and the fields E\ 
and E 2 are described by Eqs. (fTS")) and (fll)]) and average 
over initial the state 1 1)3 1 0) 4 . We consider the general 
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to photodetector 



FIG. 8: (Color online) The light emitted by SPDC experiment 
is collected from a small portion of solid angle. 



case: only one of two modes /?+ and /3_ is obtained, such 
as /3+, then the peak intensity of the field Ip x is 



IeAL) 



thuja 



2e A cf [CT p 



Ui\ 2 + \A 



2| 2 lle^ L | 2 , 



(35) 



and the peak intensity of the field Ie 2 under the same 
condition is 



e A eS cT p 



2| /3+L|2 



(36) 



The ideal efficiency r/j (j — 1,2) for conversion of power 
from the pump photons tOj (j — I. 2) to the photons lo 2 +j 
(j = l,2)is 



P 3 (L) _ I El (L) 
Pi{0) /Pi(0)' 



p 2 (o)- 



(37) 



Considering the efficienc y o f the present available single- 
photon source rj s ss 9% [23j |. we can have the total con- 
version efficiencies in our scheme 77 tot . (J = 1, 2) is 



7?tot, 



VjVs, j = {l,2}. 



(38) 



Now, we numerically evaluate the total conversion 
efficiency rytotj ■ As a numeric example with a par- 
ticular 87 Rb atomic ensemble using, we assume the 
hyperfine levels involved: {|5 2 <Si/2, -F — l,mp = 
1), \5 2 S 1/2 ,F = 2,m F = 1), \5 2 P 1/2 ,F> = 2,m F , = 0), 



\5 2 P : 



3/2: 



, F' = 2, mp/ = 2)}, which correspond to the lev- 
els {|1), |3), |2), |4)} in our scheme. As a result, the 
transition dipole matrix elements \(i2i\ = 2.992 eao/'\/T2 
and I/U43I = 4.227 eao/y/12, and the excited transitions 
1 2) — > 1 3) and 14) — » |1) have the same Clebsch- Gordon 
coefficient 1 /2 [24| . The spot sizes of the probe and FWM 
laser beams u>q are assumed to the same with wq ~ 10 
/im in our calculation. The dimcnsionless two-photon de- 
tuning is chosen and At p = 3 corresponding 
to f}+ mode and the other parameters are chosen as the 
same as those using in the Fig. [3l Finally, we have the 
total conversion efficiencies for the atomic ensemble, 



r?tot a » 5.919 x 10- 8 /cm, 



and 



Vtot 2 



4.819 x 10 _9 /cm. 



(39) 



(40) 
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FIG. 9: Phase mismatching induced by Stark shift Afcshift as 
function of detuning ujt p for Ar p = — 1. The other parameters 
are chosen as those using in Fig. [3] 



For comparison, we also examine the conversion efficiency 
for the generation of photon pairs in SPDC case. A 
SPDC process is described by the Hamiltonian, H = 
ihr)A p a\a\+}i.c., with a strong classical field A p and the 
weak signal (idler) field a \ (a\). The output state in 
SPDC case can be approximated as \ip{t)} = |0) + <7|l s , 1,) 
for g — r/Apt << 1. The conversion efficiency of pump 
photons into correlated photon pairs integrated over all 
emission directions is on the order of 3 x 10 -8 mm _1 sr _1 
for a typical nonlinear material [25L [26| . However, in 
practical operation, one photodetector can only collect 
photons from a small portion of solid ang le as shown in 
Fig. [51 For example, in the experiment [2q|. the collection 
angle is 3.3 x 10 -5 sr. This leads to a realistic conver- 
sion efficiency is about 1 x 10~ 12 mm" 1 of crystal length. 
Taking account of the reduction of efficiency due to col- 
lection angle in real detection for SPDC case, it is evident 
that the overall conversion efficiency in our case will be 
higher due to the determined direction of emitted pho- 
tons by stimulated process. In this sense, our scheme can 
have some merits over SPDC case and provide a different 
method to generate photon pairs using atomic system. 



V. APPLICATIONS OF (2+1) PHOTON 
SOURCE 



In the section, we explore the possibilities for ap- 
plications of our scheme. Firstly, the scheme can be 
used to generate photon pair due to the existence of 
the second term on the right side in Eq. (|25p. This 
term represents the situation of simultaneous existence 
of two photons in the probe frequency mode (uj^) and 
one photon in the FWM frequency mode (uJi). Detec- 
tion of a trigger photon (0^3) will project this state into 
10) = aio|0) 3 |0) 4 + e-^a 2 i|l) 3 |l)4- In the weak gain 
limit considered in the paper this state \<p) has a large 
vacuum component aio and a small admixture of the two 
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photon state a^i- Such a state shares a great similarity 
to that generated in SPDC experiments. Despite such a 
similarity between our scheme and the SPDC case, our 
scheme in generating photon pairs has its own reason 
for specialty. For example, the frequency of generated 
photon pairs can be widely tuned. The conversion ef- 
ficiency can be achieved effectively higher than SPDC 
case due to the injection-seed mechanism and stimu- 
lated Raman process. The injection-seed mechanism en- 
sures a highly directional photon emission which can in- 
crease the collecting efficiency of photons in the detec- 
tion, compared to SPDC case. Stimulated Raman pro- 
cess somehow increases the emission probability. This is 
because, for an TV-photon state input, we have the term 
alaj\N) k \0)i = s/N + l\N+l) k \l)i. As a result, the pho- 
ton emission probability is N + 1 times that of the spon- 
taneous emission described by (a^,a||0)fe|0); = 
[271 ]. Though the enhancement is low for our case lim- 
ited by N — 1 , the generation probability is 2 times of 
that from spontaneous emission. Furthermore, the one- 
photon emission enhancement due to stimulated emission 
was observed by Lamas-Linares et al. [28l | . On the other 
hand, what is more important is that if a perfect sin- 
gle photon is on demand, the conversion efficiency from 
pump photons into photon pairs will be significantly im- 
proved in our scheme. Hence our scheme can have the 
potential to take over the SPDC case if the single-photon- 
on-demand source is available in the future. 

In additional to the generation of photon pairs, our 
(2+1) photon source can also be used to generate photon- 
number Fock state for quantum communication and 
other applications. Recently, Ou [27j presented a scheme 
with a number of ideas for efficient conversion between 
photons. In the scheme, one of the ideas is to use the 
four-wave mixing as a two-photon annihilator, the input 
of which requires a two-photon Fock state. Our (2+1) 
photon source can be a candidate for providing such a 
two-photon Fock state by projecting the state \ip(t)} de- 
scribed in Eq. (f2"5|) into the state |2) 3 with the detection 
of the trigger photon 1 1) 4 . 



VI. DISCUSSION 

In this section we first discuss these ac Stark-type fre- 
quency shifts. From Eqs. © and (fT0"|) , we obtain the 
phase mismatching induced by Stark shift 

Ak shlit = ReiD^Lu) - D 4 (lj)}. (41) 

Figure [9] shows phase mismatching induced by the Stark 
shift versus lot p by numerical calculations. It is clear from 
a comparison between Figs. 1^1 and [TU1 that the mismatch- 
ing Afc shlft induced by Stark shift is very small with the 
exclusion of two maximal phase mismatches that occur 
in the regions where the gain coefficients are maximized. 
From Fig. [pj we also obtain that the magnitude of two 
maximal phase mismatches are small compared with that 
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FIG. 10: Frequency-dependent gain coefficient for the probe 
field as function of detuning lut p for Ar p = —1. The other 
parameters are chosen as those using in Fig. [3] 

of the wave vector of light field, so the angle of departure 
of field E2 from the direction determine by Ak ■ z 1=3 is 
small. 

Once the single-photon probe field is introduced, we 
can generate a correlated photon pair with probability 
|«2i | 2 /2 when only one photon 1 1 ) 3 is detected. In or- 
der to obtain the desired result, the Raman gain (by the 
pump field Rabi frequency fli) should be kept very low to 
ensure the stimulated generation of one probe photon. Of 
course, the pump field Rabi frequency Q2 should be kept 
appropriately to ensure the unit conversion from state |3) 
to state 1 4) which leads to the simultaneous generation 
of one FWM photon. Figure O depicted the Raman gain 
in the Fourier transform space for the parameters used. 
As expected, the gain is small and rather flat as required 
under the condition that the ut p is small. 

We further note that with sufficient O2 our scheme has 
approximately a EIT-likc behavior, as can be seen from 
Hamiltonian ([TJ. In this case, our four-level system is 
a hybrid system which is composed of an effective EIT 
(for the FWM photon generation) and Raman gain (for 
the probe photon generation). Specifically, consider the 
limiting case where £4 = 0, A = 0, = |f2i| = l^l, 
and very large one-photon detuning 8\, we find that the 
group velocity 1/V g = 1/V g+ = 1/V g _ w l/c+K/(2\n\ 2 ) 
which is EIT-like. 

Next, we discuss the bandwidth of the Raman gain 
and the spatial width of single photons wave packet. 
From Fig. (TUJ we know the magnitude of Raman gain 
bandwidth is 10 MHz. The injecting seed single-photon 
wave-packet should be microsecond or sub-microsecond 
pulse. Such single photon source has so far only been 
achieved with radiating object (atoms, quantum dots) in 
high-finesse microcavities because the efficiency of sin- 
gle photon source is hard to obtain high in free space 
due to the light-collecting lens covers only a fraction of 
the full 47T solid angle. Using a single atom makes it 
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possible to produce single photons with controlled wave- 
form [23|, [2t| [33, [3l[ and polarization (32j, which al- 
lows realizing deterministic protocols in quantum infor- 
mation science (33[. The deterministic and high efficient 
single-photon source [H, [HI, [H, [H[ can support the 
present scheme. The single-phot on-g eneration probabil- 
ity is about 9% by Rempe group 23]. The low efficiency 
of single photon sources degrades our scheme. When a 
perfect single source is on demand our scheme will have 
potential applications. 

Finally, we point out that the effect of quantum noises 
from the atomic system is ignored since the terms of 
Langevin-noise operators are excluded in Eq. ([2]) . The ig- 
norance of quantum noises is valid here as we exploit only 
the low-gain regime of the active-Raman-gain medium for 
the purpose of generation of correlated photons. 

VII. CONCLUSION 

In conclusion, we have studied a (2+l)-photon genera- 
tion scheme using a life-time broadened four-state atomic 
system. This scheme is based on a double-A excitation 
configuration with the first A branch forming an active- 
Raman-gain medium. This is very different from the 
conventional SPDC scheme. Furthermore, it is different 
from the spontaneous emission based biphoton genera- 
tion scheme because of the injection-seeding mechanism. 
This injection-seeding technique leads to the highly di- 
rectional generation of desired photons, resulting in high 
detection efficiency. An important feature of the present 
scheme is that two identical probe photons, because of 
the stimulated Raman emission process due to injection- 



seeding, and a FWM photon are generated simultane- 
ously, yielding correlated and entangled (2+1) photons. 
Consequently, one of the probe photons can be used as 
a coincidence trigger whereas the remaining probe and 
FWM photons form a correlated or entangled pair that 
allows further experimental studies of paired propaga- 
tion. In addition, (2+l)-photon source can also be used 
to generate photon-number Fock state. Hence if the 
single-photon-on-demand source is available in the fu- 
ture, our scheme can have the potential to take over the 
SPDC case and have potential applications. 
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